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Fourier-transfrom Ghost Imaging for pure phase object based on Compressive 

Sampling algorithm 

Hui Wan^Q and Shensheng Har{j] 
Key Laboratory for Quantum Optics and Center for Cold Atom Physics of CAS, 
Shanghai Institution of Optics and Fine Mechanics, 
Chinese Academy of Sciences Shanghai 201800, China. 
(Dated: November 20, 2009) 

A special algorithm for the Fourier-transform Ghost Imaging (GI) scheme is discussed based on 
the Compressive Sampling (CS) theory. Though developed mostly in real space, CS algorithm could 
also be used for the Fourier spectrum reconstruction of pure phase object by setting a proper sensing 
matrix. This could find its application in diffraction imaging of X-ray, neutron and electron with 
higher efficiency and resolution. Simulation and experiment results are also presented to prove the 
feasibility. 

PACS numbers: 42.52.Ar, 42.50.Dv, 42.30.Wb,42.25.Kb 



The diffraction imaging of X-ray, neutron and elec- 
tronprovide significant methods to reveal microstructure 
0, 0, 0, H 0,11, Q. However, the traditional diffrac- 
tion imaging by X-ray, neutron and electron in ther- 
mal state could only be used in periodic structure imag- 
ing. As for the aperiodic structure at the order of the 
source wavelength, according to the traditional wave the- 
ory, the coherent radiation source, such as Free Elec- 
tron Laser (FEL), is essential for the imaging. However, 
since neutrons and electrons are fermions, their coher- 
ent sources with high brightness are in principle unavail- 
able; although it's possible to obtain a coherent sources 
of photons with high brightness, the high requirement for 
brightness could only be met on the Synchrotron Radia- 
tion Facilities. 

It has been proved that two-photon correlation Ghost 
Imaging (GI) could realize the diffraction imaging by tak- 
ing advantages of thermal source [8[, however, it is also 
limited by the long acquisition time, as the correlation 
theory calls for mass samples to guarantee the ensemble 
average. Actually, lots of methods were introduced to 
improve its Convergence [9j], however, as long as it still 
takes the correlation algorithm, the concept of ensemble 
makes it less effective. 

Recently the Compressive Sampling (CS) algorithm at- 
tracts more and more attention because of its extraordi- 
nary effect of reducing the samples [TO, EH EIL EE EH EH • 
GI and CS have similarities in extracting imaging: both 
are prepared with random sensing signal to "express" 
the imaging and a "point" detector to collect the result 
of expression. But they also have intrinsic difference: 
GI is based on accurate "point by point" measure model 
[H, E3> while CS theory has proved that "global ran- 
dom" measure model exhibits hig her efficiency in imaging 
extraction 0, EH, El El, El, EB| . Therefore, by combing 
GI and CS, it's possible to develop a brand-new imaging 



model with higher imaging efficiency and resolution. 

At present, CS theory concentrates mostly in the real 
space, which limits its app lication in the treatment of 
grayscale image [To|, EH, E2, E3, EI El- While the actual 
information media waves contain not only amplitude but 
also phase information, thus the existing CS algorithm 
could not be used directly to arbitrary GI scheme, es- 
pecially in the Fourier GI scheme, where most spectrum 
information relies on the phase part of the field. 

In view of these problems, we firstly propose a combi- 
nation of the Fourier GI and CS in this letter. A recov- 
ery algorithm is investigated for the Fourier- Transform 
Ghost Imaging with Compressive Sampling (GICS) with 
simulation and experiment results. This algorithm could 
not only improve the imaging efficiency to reduce the ra- 
diation damage of the sample, but may also improve the 
spectrum resolution compared to the Correlation Ghost 
Imaging (CGI). Therefore GICS may provide a brand- 
new diffraction microscopic imaging technic of X-ray, 
neutron and electron for the aperiodic structure. 

According to the CS theory, the algorithm is based on 
a corresponding relation between the imaging informa- 
tion X and the detect si gnal Y through a proper sensing 
matrix A 0, EH, E2, E1J3 : 



AX = Y 



(1) 



If we sample K times, A should beaKxiV matrix, and 
Y should be a known K-element vector. Then by solving 
a convex optimization program, it gives an optimal result 
of X with N-element in the same expressing space as A: 



min \\xWft subject to Y = AX 



(2) 
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Where x are the elements of X in certain express basis, 
and \\x\\ h = EM- 

In our case, the lensless Fourier-transform GI scheme is 
shown in Fig 1 [8] . The thermal field from the source S is 
divided by the beam splitter (BS) into test and reference 
arms. To perform the Fourier-transform GI, there should 
be d\ = c?2i + for the scheme. Based on this scheme, 
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FIG. 1: The lensless Fourier-transform GI scheme. Beam 
splitter (BS) is introduced to divide the field into two arms: 
in the test arm, field propagates freely cfoi to an object and 
then c/22 to array detector D2] in the reference arm, field 
propagates dl to another array detector D\.d\ — c/21 + ^22- 



we could express the intensity I w (r 2 ) on the test detector 
D2 as: 



I w (r 2 ) oc J dxdx'E*(x)E(x')t*(x)t(x') (3) 

obj 

x exp j-^- [(x - r 2 ) 2 - (x' - r 2 ) 2 ] | 

Where the integration is over the object plane. By 
comparing with the CS theory, the object information is 
located on the right side of ([3|), thus I w {r 2 ) could be con- 
sidered as the known detect signal Y. Since our imaging 
goal is the Fourier-transform of t(x), to process the CS 
algorithm, it's necessary to establish a corresponding re- 
lation like (pQ), and the sensing matrix A must be related 
to the signal I r (r\) from the reference detector D\ as the 
spectrum used to be expected to show up on the array 
detector D\{r\). 

If we consider the array of D\ to be large enough, then 
(|3|) could be rewritten approximately as 



,(r 2 ) oc J dndri£7*(n)£7(ri) j dxdx' exp {-^- [(x - n) 2 - (x' - r[) 2 ] J 
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(4) 
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(5) 



Compare (j4j) with ([T]), the sensing matrix A could then 
be expressed as E*(ri)E(r[) exp |-^ 7r /A(i 22 ( r i _ r i 2 ) }• 

Worth to note, unlike the sensing matrix A introduced 
in the bucket real-space GI scheme [15] . the matrix A 
here is a K x N 2 one rather than a K x TV one. 

This arouses two problems: one is that the CS theory 
is normally developed in the real space, while the rela- 
tion (|4]) established in the complex space is not directly 
available in CS algorithm; the other is the unobtainable 
interference terms of A when n 7^ r[. 

To solve the first problem, A should be redesigned. 
Define A = Al+iA2 and the unknown spectrum T*((ri — 
r 2 )/(Ad 22 ))r((ri - r 2 )/(Ad 22 )) = B x + zB 2 . Obviously, 
Al, A2 and 51, 52 are all real matrixes with dimensions 
of if x N 2 and iV 2 x 1 respectively. Then (|4j) could be 
replaced by two equations: 



A X B X - A 2 B 2 = I w (6) 



A X B 2 +A 2 B X =0 (7) 

It seems there're IN 2 degree of freedom for X until we 
realize the symmetry of Al and the anti-symmetry of A2. 
This in return guarantees the symmetry of Bl and the 
anti-symmetry of B2 according to the convex optimiza- 
tion ([2]) in CS progressing, and make ([7j) always hold. 
Therefore we can rewrite the new sensing matrix A 1 by 
combing the elements of Al and A2 without introducing 
extra degree of freedom beyond TV 2 . 

Fig 2 is an example to rebuilt A' for a single sampling: 
its upper triangular elements are the same as those of 
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FIG. 2: An example to rebuilt A' . 



2 x Al, while the lower triangular ones are the same as 
—2 x Al, and the diagonal line ones are the square root 
of N intensities y/I r (r\) from N points on Dl. Thus by 
performing the CS algorithm, as long as we know the 
exact field distribution on Dl, not only the spectrum in- 
tensity |T (/)| is obtainable on the diagonal line, but also 
its phase information from the non- diagonal elements. 

Generally, the only way available to get the exact field 
distribution on Dl in GI scheme is by calculating for 
the controllable source scheme or by introducing another 
known reference field E' for the homodyne detection. 
However, both are unavailable for the diffraction imaging 
of X-ray, neutron and electron. Notice the final imaging 
spectrum result obtained from correlation GI scheme is 
just |T(/)|, which lies on the diagonal line of the CS 
solution X, corresponding to the known intensity distri- 
bution of Dl on the diagonal line of A! . According to the 
CS algorithm, the relative distribution of |T (/)| could be 
expressed through the known l r {r\) of K sampling, while 
the existence of other non-diagonal elements just make 
the relation (j6]) hold. Therefore, if we only want to get 
|T(/)|, the accuracy of the non-diagonal elements' val- 
ues, mainly influencing the convergence, become less im- 
portant. And it's possible for us to proper conjecture the 
phase distribution on Dl based on the concrete scheme 
and translate it into A' followed by normal CS algorithm. 

Obviously, (|6|) is not a strict equation in the phase con- 
jecture course, just as a non-ideal bucket detector for the 
real-space imaging [l5|, whose influence to the sensing 
matrix is ignored. However, this approximation only in- 
fluences the efficiency to extract the information rather 
than destroying the expressing basis. Also, there's little 
relation in applications with strict equation for the CS al- 
gorithm to perform, thus it makes little sense to pursue 
the strict condition for CS algorithm. On the contrary, 
it's proper to develop the CS algorithm adaptive to rela- 
tions with different approximation. And still, GICS, as 
a brand new imaging algorithm instead of CGI, brings 
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FIG. 3: Transmission function and spectrum intensity \T (/)| 
of the pure phase five-slit object, (a) Transmission function of 
the object, the slit wide a is 600 /mi. (b) Spectrum intensity 
|T(/)| of the object observed from a lens of 30 mm focal 
length. 



improvement of the low efficiency and limit resolution to 
the traditional CGI. 

Our simulation scheme is the same as shown in Fig.l, 
with source size to be 3 mm, cfei = 20cm, ^22 = 
20cm, di = c?2i + ^22, both Dl and D2 are array de- 
tectors. The object is a one dimension pure phase five- 
slit, whose transmission function and spectrum intensity 
\T (f) \ is shown in Fig 3. We simulate the result of GICS 
with homodyne detection and phase conjecture courses 
in Fig 4 (a) and (b). For comparison, we also simu- 
late a result of GICS algorithm applied merely to the 
square root of N intensities from Dl (shown in Fig 4 
(c)) to show the improved convergency by introducing 
in the non-diagonal elements. All results are performed 
through spatial averaging treatment over K=50 sampling 
Q. We could see that GICS could also extract the de- 
sired information with higher efficiency even without a 
strict equation (j4]). Besides, the extra non-diagonal ele- 
ments, though just coming from conjecture, improve the 
extract efficiency. 

We also take experiments of the same scheme, where 
a pseudo-thermal light is produced by permeate a laser 
beam with diameter of 3 mm through a rotating ground 
glass, c?2i = 20cm, ^22 = 5cm, d\ = di\ +^22, and the ob- 
ject is a similar one dimension pure phase five-slit with 
a = 150 jam. The sampling number K = 100, and the re- 
sult has been performed through spatial averaging tech- 
nique. 

In conclusion, A new algorithm based on the lensless 
Fourier- transform GICS is presented in this letter. This 
scheme opens up a new way for the diffraction imag- 
ing of X-ray, neutron and electron with higher efficiency 
and resolution. The GICS, however, is much more than 
a algorithm revolution. From the point of microscopic 
essence, quantum fluctuations of field are generally de- 
scribed via high order correlation, which is based on ac- 
curate classical "point by point" measurement: the tra- 
ditional imaging technic is first order correlation (single- 
photon) system, while CGI is a second order correlation 
(two-photon) system. CS, on the contrary, stems from a 
"global random" sensing measurement. By introducing 
CS algorithm into GI, it also initiates a new discussion 
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FIG. 4: Simulation results of GICS with homodyne detection, phase conjecture courses based on the Fig 1 scheme with 
d>2i =, (I22 =, c?i = cfoi + o?22- (a)result of GICS with homodyne detection .(b)result of phase conjecture courses. 




FIG. 5: Experiment result of GICS based on the GI scheme 
shown in Fig 1, where di\ — 20cm, dii — 5cm. The 



on its microscopic quantum mechanism. To find a rea- 
sonable explanation for the microscopic essence of GICS 
is a challenging problem, which will help to develop new 
imaging model with higher efficiency. 
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